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FORECASTING  THE  FOCUS  OP  AIR  BIASTS  DUE  TO ' 

METEOROLOGICAL  CONDITIONS  IN  THE  LOWER  ATMOSPHERE 

ABSTRACT 

Whenever  explosions  are  used  in  testing  or  in  experimental  ’ 
procedures,  the  sound  waves  that  go  beyond  the  limits  of  the  installation 
may  cause  ccxaplaints  of  annoyance  or  damage  from  otherwise  good  neighbors 
This  is  due  to  focusing  of  the  soxind  waves  caused  by  the  meteorological 
conditions  at  the  time  of  the  explosion. 

The  theory  of  the  propagation  of  sound  through  the  atmosphere 
is  given  briefly.  The  conditions  in  the  atmosphere  which  cause  the 
sound  to  be  focused  are  velocity  gradients  produced  by  variations  with 
altitude  of  humidity,  air  temperature,  and  wind  velocity.  A  simple 
method  is  described  for  evaluating  these  factors  and  forecasting  the 
location  of  a  focus,  if  one  is  to  be  expected,  as  well  as  the  intensity 
of  the  sound  at  the  focus. 

This  technique  has  been  employed  successfully  at  Aberdeen  Proving 
Ground  for  three  years  during  which  time  the  complaints  have  been  greatly 
reduced. 
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BACKGROUND 


At  every  Department  of  Defense  establishment  where  testing 
operations  Involve  the  explosion  of  H.E.  charges  outdoors^  there  are  . 
frequent  complaints  of  annoyajice  or  claims  for  damage  from  persons 
living  in  contiguous  areas.  The  complaints  of  annoyance  are  often 
justified;  the  damage  claims  are  often  unjustified.  Both  can  be 
reduced  or  avoided  entirely. 

The  refraction  of  soiuid  waves  produced  by  the  meteorological 
conditions  between  the  earth's  surface  and  the  elevation  of  10^000  to’ 
12^000  feet  may  cause  the  somd  waves  produced  by  explosions  at  or  above 
the  surface  of  the  ground  to  be  focused  near  homes  in  the  vicinity  of 
the  testing  area.  This  phenomenon  may  be  annoying;  however,  it  is 
doubtful  that  any  damage  to  walls  or  window  glass  will  be  caused  by  air 
blast  at  distances  greater  than  a  few  miles  from  any  but  very  large 
explosions  (100  or  more  pounds  of  high  explosive). 

.  The  refraction  of  the  sound  waves  is  due  to  the  presence  of  vertical 
gradients  of  the  velocity  of  sound  through  the  air  (^)  where  v  is  the 
velocity  of  sound  and  y  is  the  altitude.  Since  the  velocity  of  sound 
in  air  depends  on  temperature,  humidity,  and  the  wind,  the  vertical 
gradient  will  depend  on  variations  of  these  factors  with  altitude. 

One  of  the  earliest  papers  on  the  refraction  of  sound  due  to 
meteorological  conditions  was  presented  to  the  Royal  Society  of  London 

•  t 

about  1906  as  a  result  of  the  firing  of  guns  during  the  funeral  of  Queen 
Victoria.  The  easily  recognized  noise  resulting  from  the  rythmic  firing 
of  the  guns  was  heard  in  London  and  also  far  to  the  north  but  the  two 
areas  were  separated  by  a  zone  of  silence.  This  experience  was  explained 
by  refraction  due  to  a  wind,  directed  to  the  north,  with  a  large  wind 
gradient  existing  at  high  altitudes.  Since  that  time,  about  sixty  papers 
have  been  published  on  this  subject.  Everett  P.  Cox^^^,  in  his  treatment 
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Nimerals  refer  to  references  listed  at  end  of  the  report. 
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of  upper  air  temperatures  as  measured  at  the  demolition  of  Helgoland,  . 
has  reviewed  the  salient  facts  and  presented  an  excellent  ‘bibliography 
of  the  important  contributions.  ; 

I 

■  In  the  present  paper,  nothing  has  been  added  to, the  classical, 
theory  of  propagation;  however,  an  operation  procedure  is  described, 
that  permits  the  meteorological  conditions  to  be  quickly  evaluated  and 
the  cause  of  complaints  from  contiguous  areas  reduced  to  a  minimimi. 


BASIC  THEORY 


As  stated  prevloualy,  the  velocity  of  sound  In  the  atmosphere  depends 
on  temperatvire ,  humidity,  and  vlnd.  E3q)erlence  shows  that  the  ten^jerature 
of  the  lower  atmosphere  at  Aberdeen  Proving  Ground  may  vary  through  a 
range  of  20  to  50°C  between  the  surface  and  an  altitude  of  10,000  feet, 
causing  a  change  In  velocity  of  sound  of  40  to  60  feet  per  second.  3n 
this  same  region,  the  wind  may  vary  from  0  to  4?  or  60  feet  per  second. 
Although  the  change  In  relative  humidity  can  be  as  large  as  50^,  such  a 
change  can  affect  the  velocity  of  sound  by  only  about  1  foot  per  second. 
Since  1  foot  per  second  is  less  than  the  error  In  our  determination  of 
wind  velocities,  the  effect  of  changes  In  humidity  can  be  neglected  without 
noticeable  error. 

When  the  problem  of  devising  a  simple  but  safe  set  of  guidance  rules 
for  firing  at  Aberdeen  Proving  Ground  was  given  to  the  Ballistic  Research 
Laboratories,  a  study  was  made  of  the  meteorological  conditions  prevailing 
in  the  vicinity  of  Aberdeen  on  those  days  when  complaints  had  been 
registered.  Of  twenty-six  cases,  in  all  but  one,  there  was  a  strong 
wind  increasing  in  velocity  at  altitudes  of  about  two  or  three  thousand 
feet  up  to  eight  or  ten  thousand  feet.  The  wind  was  in  general  from 
the  point  of  detonation  toward  the  point  of  complaint.  This  fact  was 
taken  as  svifficlent  evidence  that  any  guidance  rules  for  the  Aberdeen 
Proving  Ground  area  should  include  consideration  of  wind  effects  in 
addition  to  the  effect  of  temperature  on  the  propagation  of  the  sound 
waves.  Thus  velocity  gradients  are  based  on  an  algebraic  sum  of  the 
temperature  effect  and  the  wind  effect  at  each  altitude .  The  temperature 
effect  Is  of  course  a  scalar  quantity  and  Is  the  same  In  all  directions, 
whereas  the  wind  effect  Is  a  vector  quantity  and  Is  a  function  of  direc¬ 
tion.  Before  describing  how  the  temperature  and  wind  effects  con  he 
calculated,  the  nature  of  the  paths  of  sound  waves  will  be  reviewed 
briefly . 

When  an  explosion  occurs  at  the  surface  of  the  ground,  a  diverging 
shock  front  starts  at  the  point  of  e^qtloslon  and  spreads  in  all  direc¬ 
tions  in  the  hemisphere  above  the  surface  of  the  ground.  The  shock 
waves  degenerate  to  sound  waves  quite  rapidly  and  it  is  sufficient  for 
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the  purpose  of  this  study  to  consider  only  the  propagation  of  sound 
waves . 


The  paths  of  the  various  parts  of  the  wave  front  can  be  depicted 
by  rays  emanating  In  all  directions  from  the  center  of  the  explosion. 
One  of  these  rays  will  start  from  the  sororce  making  an  angle  6  with  the 
horizontal.  Then  at  the  surface 


y’ 


dx 


tan  0 


where  x  and  y  are  the  horizontal  and  vertlcle  coordinates.  If  0 

y 

the  angle  betveen  the  ray  and  the  horizontal  at  any  eiltitude  y; 


is 


and 


y*  = 


^  =  tan  0 
dx  3 


y  = 


2  d  (tan  0  ) 
d  y  '  y' 

dx 


(1) 

(2) 


If  the  velocity  of  sound  varies  with  altitude,  the  angle  0  will  vary 
according  to  Snell’s  law  such  that.  If  V  is  the  velocity  at  altltxode 

y 

y,  then 

coa.0  . 

— -5 — =  a  constant  =  C  (5) 

y 


or 


from  (l) 


from  (4)  and  (5) 


or 


Cos  0y  =  CVy 


sin  0 


(M 

(5) 

(6) 

(7) 
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and 


-  dV 

V’  -  -  y 

^  '  2  3  dy 

C  V 

y 


y" - -  ^ 


dV 


^  dv 
V  cos  0  ^ 

y  y 


(sec^  9  ) 

Vy  ^  y'  dy 


dV 


(8) 

(9) 


=  ^[i-(y)]^  (10) 

'  dv 

If  =  a  constant  K  (the  vertical  velocity  gradient)  and  If  »  velocity 
at  the  ground  surface,  then 

Vy  =  V^+’Ky^  (11) 


=  K  ( 


+  y  ) 


(12) 


and  y"  =  - 


The  solution  of  equatipn  13  (see  references  1,  2,  and  5)  represents  the 
path  of  propagation  with  least  time  "between  two  points:  ' 

(x  +  C^f  +  (y  +  (14) 

In  which  and  Cg  are  constants  of  integration,  and  x  and  y  are  the 
cartesian  coordinates  of  the  points  along  the  ray  i)ath.  Equation  l4 
represents  a  circle  the  center  of  which  is  on  a  line  helow  the  groxmd 
surface  a  distance  eqi^l  to  V^^/K.  The  parameter  Is  the  radius  of 


K  (  ^  +  y  ) 
K  _ 


1  +  (y*)" 

V 


1  +  irf 


K 


1  +  y 


K 


(13) 
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the  circle  and  is  the  horizontal  distance  fraii  the  origin  of  co¬ 
ordinates  to  the  center  of  the  circle.  In  Figure  1  are  shown  the  paths 
of  several  rays  starting  at  the  source  with  different  angles  of  departure. 
For  a  negative  gradient  the  path  of  the  sound  wave  would  be  the  arc  of  a 
circle  curving  toward,  the  center  of  the  circle  being  on  a  line  above  the- 
ground  surface  at  a  distance  equal  to  V^^/K. 

Variations  with  altitude  of  either  the  rate  of  change  in  air  tem¬ 
perature  or  the  rate  of  change  in  wind  velocity  will  change  the  value  of 
the  gradient  K  and  therefore  the  curvature  and  direction  of  the  cor¬ 
responding  arc  will  change  accordingly.  At  the  boundary  between  layers 
having  different  values  of  K  there  will  be  a  layer  of  transition  in  which 
the  conditions  gradmlly  change  from  those  of  one  layer  to  those  of  the 
other.  In  this  transition  layer  the  path  of  a  ray  gradually  changes  from 
the  arc  determined  by  the  velocity  gradient  in  the  layer  from  which  the 
ray  is  emerging  to. the  arc  determined  by  the  gradient  of  the  layer  being 
entered.  Figure  2  shows  a  typical  graph  of  sound  velocity  versus  alti¬ 
tude  and  Figure  3  shows  the  paths  of  the  rays  resulting  frcm  the  gradients 
in  Figure  2.  As  indicated  above,  the  curvature  of  the  paths  in  any  layer 
depends  on  the ' gradient  in  the  layer . 

For  a  ray  to  be  refracted  to  the  surface  from  any  layer,  the  maximum 
velocity  attained  in  that  layer  must  equal  or  exceed  the  velocity  at  all 
points  in  the  medium  nearer  the  surface  of  the  ground. 

If  the  ray  path  is  confined  to  a  medium  in  which  only  a  single 
velocity  gradient  exists  the  horizontal,  distance  traversed  by  a  given 
ray  (l.e.  the  range)  can  be  most  simply  calculated  by  plane  gecanetry. 

Ih  Figure  4  a  ray  path  is  shown  leaving  the  origin  at  an  angle  0,  rising 
in  an  arc  to  an  altitude  y  and  returning  to  the  surface  at  a  horizontal 
distance  R  from  the  origin. 

In  the  figure. 


FIG,  I  -  PLOT  OF  EQUATION 
i.€.  FOR  DIFFERENT 
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LtNE  OF  CENTERS 
FOR  POSITIVE  GRADIENT 


14  FOR  VARIOUS  VALUES  OF  C| 
ANGLES  OF  DEPARTURE 
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ALTITUDE  IN  KILO- FEET 


FIG  2-TOTAL  CHANGE  IN  VELOCITY  vs  ALTITUDE 
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1 


ALTITUDE  IN  KILO  FEET 


FIG.  3- RAY  PATHS  CORRESPONDING  TO  GRADIENTS 
REPRESENTED  IN  FIG.  2 
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B  =  ANGLE  OF  DEPARTURE 
R  -  HORIZONTAL  RANGE 

y  =  MAXIMUM  ALTITUDE  REACHED  BY  THE  RAY. 


FIGURE  4.“  RANGE  AND  ALTITUDE  FOR  A  GIVEN  RAY  IN 
A  MEDIUM  OF  A  SINGLE  GRADIENT 
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where  =  initial  velocity  and  K  is  the  velocity  gradient. 


The  travel  time  t  =  2 


Jq  Y 


In  which 
and 

After  integration 


ds  =  +  (iy)^ 


1/2 


Vy=V^+Ky 


t  =  I  log  cot  (  -  I  ) 


K  6 


For  a  medim  having  two  gradients  the  path  of  a  ray  is  pictured 
In  Figure  5.  The  range  of  the  ray  departing  at  an  angle  0^  will  be; 

R  =  2X^  +  2  Xg 
=  a  c  -  b  c 

ac  =  ^  tan  0, 


tan  0^ 


V2  V, 


Since  Vg  -  \  - 


\  ^2 

and  ■tan  “tan  0g 


Xg  .  ^  tan  9g 


2V, 


and  R  = 


.2 

^2 


2V  sin  9 

tan  0^  +  — i  - - 

2  cos  0j^ 


^Vg  sin  ©2 


cos  9,. 
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=  ANGLE  BETWEEN  RAY  AND  HORIZONTAL  AT  SURFACE. 

$2  =  angle  BETWEEN  RAY  AND  HORIZONTAL  AT  TOP  OF  LAYER, 
a  d  =  RANGE  OF  RAY=  2X,  + 

Kj.Kg  =  GRADIENTS  IN  Ist  AND  ZQd  LAYER  RESPECTIVELY. 


FIGURE  5  -  RANGE  OF  RAY  IN  A  MEDIUM  OF  TWO  LAYERS 
HAVING  VELOCITY  GRADIENTS  OF  K,  a  Kg 
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vhich  jreduces  to- 


2Vp  2  V, 

E  -5^  tan  62  ^  (=ln  6^  -  sin  6^  )  (20) 

The  trasrel  time  iCor  the  two  layered  medim  Is 


/y  ^  y 

^  ^  (for  1st  layer)  +2  j  ^  ^  (for  2nd  layer) 

n  ’v  •/v  '^v 


For  1st  layer 


y 


ds  r  \  2  2 

^  in  which  ds  =  (dx)  +  (dy)  and  V  =  +  Ky 

''y  L.  — I  y  . 


/n 


t  -  2  f  ^  (  sin  9  )  dy  _  2 

"1-  I  V,  -.Ky  -  7 


Qr 


-  do 


cos  9 


9 

log  tan  (  ^  +  ^  ) 


H- 

«1 


r 


©1 

n  .  1 


J»1 


Ic 

_ 


©■1 

n  1 


log  cot  (  ^  -  ^  )  -  los  cot  ( 


i‘r  ’] 

>] 


«  ®2 
¥  ■  2“ 


For  2nd  layer 


t2=2 


•/  .  V 


in  which  +  K^Y 


=  ^  log  cot  ( 


ff  ®2 
¥  "  ?" 


The  total  travel  time  t  =  tj^  +  t^ 
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(21) 


t 


|-  |log  eot  <  J  )  - 

+  ^  log  cot  (  5  ~  2^  ^ 


log  cot  ( 


The  symbols  of  equation  (2l)  are  defined  in  the  diagram  of  Flgin?e  6. 

Figure  7  presents  a  plot  of  the  travel  time  (t)  versus  angle  of 
departure  (q)  and  range  (r)  versus  Q  for  the  case  of  transmission  throvigh 
an  atmosphere  in  vhich  two  velocity  gradients  and  exist,  being 
less  than  K^.'  It  should  be  noted  that  at  0  =  the  travel  time  curve 
and  the  range  pass  thru  a  minimum.  Furthermore,  both  the  travel  time 
and  range  change  very  slowly  with  ©  in  the  vicinity  of  ©j^.  This  indicates 
that  for  the  conditions  stated  for  Figure  7)  s-  cone  of  rays  will  be 
converged  at  arrange  of  with  very  slight  differences  in  phase  which 
will  produce  an  Increase  in  Intensity.  The  various  combinations  of 
gradients  that  will  produce  this  focussing  is  of  vital  Interest  and  can 
best  be  shown  by  calculating  ray  paths  for  various  combinations  of 
gradients  that  can  be  expected  to  occur.  An  analogue  computer  has  been 
used  to  calculate  and  trace  ray  paths  for  single  gradients  and  for 
combinations  of  two  to  seven  gradients.  The  solutions  show  paths  for  a 
group  of  rays  starting  at  different  angles  of  departure.  If  focussing 
occurs  the  convergence  will  be  shown  as  well  as  the  distance  to  the 

f 

point  of  focus. 


The  Piedmont  Division  of  ^erry-Rand  Corporation,  CharlottesvliJle, 
Virginia, designed  the  analogue  computer  used.  With  it  the  effect  of 
from  one  to  twenty  gradients  can  be  considered.  The  computer  is  small 
in  volume  (about  two  cubic  feet)  and  provides,  the  solution  in  a  few 
minutes.  In  order  to  simplify  the  construction  of  the  computer  some 
simplifying  assumptions  were  made  in  the  solution  of  the  equation  of  the 
ray  path.  The  approximations  and  the  resulting  errors  are  Insignificant. 
This  is  demonstrated  in  Appendix  A. 
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PATH  OF  RAY  THROUGH  TWO 
LAYER  ATMOSPHERE 


FIG  6 -DIAGRAM  DEFINING  TERMS  IN  EQUATION  21 
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BLAST  WAVE  HJTENSm 


The  overpressure  In  a  blast  wave  when  propagated  through  an 
atmosphere  In  which  the  velocity  gradient  is  zero  has  been  determined  , 
for  various  distances  for  several  different  weights  of  explosive.  At 
a  given  distance  the  overpressure  in  the  blast  wave  will  be  different  : 
if  the  explosive  is  detonated  at  the  surface  of  the  ground  or  in  free 
air.  The  measured  value  will  also  be  different  if  the  pressure  gage  is 
at  the  surface  of  the  ground  or  In  free  air.  In  Figure  8  graphs  are 
shown  which  give  the  overpressure  versus  distance  for  several  weights  of 
explosive  detonated  at  the  surface  and  the  overpressure  measured  at  the 
surface  for  a  uniform  atmosphere. 

When  the  ray  paths  of  sound  are  affected  by  gradients  in  the 
atmosphere  the  Intensity  at  a  given  distance  will  be  some  multiple  of 
the. overpressure  given  in  Figure  8  for  the  given  distance.  The  com¬ 
binations  of  gradients  to  be  experienced  can  be  classed  as  one  of  five 
basic  types  or  categories.  The  types  are  shown  in  Figure  9*  For  each 
type  approximate  values  have  been  determined  for  a  ''multiplication 
factor"  to  calculate  the  intensity  at  a  focus  due  to  that  type  of 
gradient  combination. 

The  value  of  intensity  at  a  focus  compared  to  the  intensity  at 
the  same  distance  in  a  uniform  atmosphere  has  been  derived  from  the 
experience  of  several  years.  A  few  direct  and  many  indirect  deter¬ 
minations  have  been  made.  By  noting  the  distance  to  a  peurtlcular  type 
of  damage  and  assuming  the  minimum  overpressure  known  to  produce  such  • 
damage,  a  maximum  "multiplication  factor"  can  be  calculated.  Such 
indirect  methods  have  provided  conservative  but  very  useful  values  for 
calcvilating  overpressures  to  be  expected.  The  calculation  of  the 
Intensity  to  be  expected  at  the  focus  Is  made  as  follows; 

The  distance  from  the  explosion  to  the  focus  is  noted  in  the 
solution  on  the  computer,  eind  the  intensity  of  the  blast  wave  at  that 
distance  in  a  uniform  atmosphere  is  read  from  Figure  8.  The  category 
of  the  gradient  combinations  in  the  atmosphere  is  determined  by  compar¬ 
ison  of  the  velocity  versus  altltiide  curve  with  Figure  9  aud  the 
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FIG.  9 


DESCRIPTION 


MULTIPLICATION 

FACTOR 


SINGLE  NEGATIVE  GRADIENT 


SINGLE  POSITIVE  GRADIENT 


0 


ZERO  GRADIENT  NEAR  SURFACE 
WITH  POSITIVE  GRADIENT  ABOVE 


10 


WEAK  POSITIVE  GRADIENT  NEAR 
SURFACE  WITH  STRONG  POSITIVE 
GRADIENT  ABOVE 
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NEGATIVE  GRADIENT  NEAR  SURFACE 
WITH  STRONG  POSITIVE 
GRADIENT  ABOVE 


100 


-VARIOUS  TYPES  OF.  VELOCITY  GRADIENTS  TO 
BE  EXPECTED  AND  THE  INCREASE  IN  INTENSITY 
AT  A  FOCUS  FOR  EACH  TYPE. 
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"multiplication  factor"  is  observed.  The  Intensity  read  from  Figure  ,8 
multiplied  by  the  "multiplication  factor"  will  be  the  intensity  at  the 
focus  to  be  expected. 

DAMAGE  CAUSED  BY  BLAST  WAVES 

The  damage  caused  by  the  overpressure  of  a  blast  wave  depends  to 
a  very  large  extent  on  the  type  of  construction.  Glass  panes  of  average 
size  and  thiclmess  vary  greatly  in  their  ability  to  withstand  blasts 
depending  on  how  the  panes  are  mounted.  If  the  pane  is  forced  into  the 
frame  so  as  to  be  under  a  constant  strain,  a  blast  wave  of  0.1  psl 
overpressure  can  cause  the  pane  to  crack  but  if  it  is  mounted  without 
any  strain  In  the  glass  an  overpressure  of  about  0.75  Psi  may  be 
required  to  crack  it. 

The  cracking  of  plaster  on  a  wall  depends  on  the  flexibility  of 
the  wall.  A  plastered  surface  attached  to  a  masonry  wall  will  withstand 
a  much  higher  pressure  than  a  surface  supported  by  a  wide  wooden  panel. 
In  general  a  well  constructed  plastered  wall  will  stand  higher  over-  ■ 
pressures  than  average  window  panes. 

An  overpressure  of  .03  to  .05  psi  in  a  blast  wave  can  cause  a 
loose  window  sash  to  slap  the  window  frame  and  produce  a  loud  noise 
while  actually  no  damage  is  being  done.  In  contrast  to  this,  the  quiet 
settling  of.  one  corner  of  a  house  can  cause  damage  to  walls  and  windows 
which  is  often  attributed  to  blast  waves. 

GATBERING  METEOROLOGICAL  DATA  •  -  . 

The  meteorological  data  are  gathered  by  release  of  a  weather 
balloon  and  radiosonde-.  ”  A  GMD  1  A  Rawin  set  is  used  at  Aberdeen  Proving 
Ground.  The  wind  data  are  calculated  from  the  successive  positions  of 
the  balloon  which  is  determined  at  frequent  intervals.  The  air  tern- 
peratxure  is  telemetered  back  by  the  radiosonde.  From  the  ground  surface 
to  5000  feet  altitude,  the  temperature  and  wind  velocity  and  direction 
should  be  determined  for  5OO  feet  intervals.  From  5OOO  feet  to  about 
12000  feet  the  data  should  be  determined  every  1000  feet.  The  data 
are  recorded  as  in  the  sample  data  sheet  (Figure  lO) .  The  air  tem¬ 
perature  is  recorded  in  degrees  centigrade.  A  change  in  temperature 
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METEOROLOGICAL  DATA 
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NOTE'. 

IN  COLUMN  3-RECORD  DIRECTION  TO  WHICH  WIND  IS  BLOWING 

IN  COLUMN  4 -ANGLE  BETWEEN  WIND  DIRECTION  AND  AZIMUTH  OF  INTEREST 

In  COLUMN  6 -(WIND  VELOCITY  AT  EACH  ALTITUDE)  X  (COS  Y  AT  SAME  ALTITUDE) 

IN  COLUMN  9- (VALUE  IN  COLUMN  I  AT  EACH  ALTITUDE)  X  2 

COLUMN  10  IS  THE  SUM  OF  COLUMNS  8  AND  9 


FIGURE  10- SAMPLE  DATA  SHEET 
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of  the  air  of  i°C  changes  the  velocity  of  sound  almost  exactly  2  feet/ > 
second,  so  the  difference  between  the  velocity  of  sovmd  at  a  given 
temperature  and  the  velocity  at  0°C  Is  determined  by  multiplying  -Uae  i 
temperature  by  2.  Since  the  effect  of  wind  varies  with  direction,  the, 
component  of  the  wind  In  the  direction  of  Interest  must  be  calculated 
for  each  elevation. 

At  each  elevation,  the  sum  of  the  component  of  the  wind  In  the 
direction  of  interest  (colxmin  8  In  Figure  10)  and  the  change  in  velocity 
due  to  the  temperature  (column  9)  will  give  the  difference  between  the 
velocity  at  that  elevation  and  the  velocity  at  zero  degrees  centigrade. 

If  the  values  of  this  resultant  (column  10)  are  plotted  on  linear  co¬ 
ordinates  as  In  Figure  11,  the  variation  in  velocity  Is  depicted  and 
the  magnitude  of  the  gradients  can  be  measured.  Figure  11,  Is  the 
presentation  of  the  data  upon  which  the  forecast  of  a  focus  and  its 
location  Is  based.  If  the  wind  varies  widely  in  direction  with  altitude, 
It  may  be  necessary  to  reduce  the  data  for  two  or  more  directions,  each 
case  being  treated  separately  as  described  above. 

FORECASTING  THE  FOCUS 

If  the  velocity  data  shown  In  Figure  11  are  furnished  to  the 
Sperry-Rand  analogue  computer,  a  presentation  of  the  resulting  ray  paths 
will  be  shown  on  the  screen.  Ray  paths  corresponding  to  the  conditions 
depicted  in  Figure  11  are  shown  In  Figure  12.  If  a  focus  occurs.  It 
will  be  seen  on  the  screen  where  the  rays  converge  and  the  distance  from 
the  source  will  be  shown. 

Actually  at  APG,  It  has  been  found  that  the  pattern  of  gradients 
often  repeat  so  that  a  library  of  32  "typical  gradients"  and  the  resulting 
ray  paths  have  been  accumulated.  Generally  the  gradients  recorded  daily 
can  be  matched  sufficiently  close  with  one  of  the  "typical  gradients" 
to  be  used  for  a  forecast.  The  computer  should  be  kept  at  hand  for  more 
precise  determination  when  needed,  and  for  those  cases  for  which  there 
Is  no  matching  typical  gradient.  The  typical  gradients  are  included 
In  Appendix  B. 
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ALTITUDE  IN  KILO  FEET 


FIG  II -PLOT  OF  TOTAL  CHANGE  JN  VELOCITY  vs.  ALTITUDE 


ALTITUDE  IN  KILO  FEET 


FIG.  12 -RAY  PATHS  CORRESPONDING  TO  GRADIENTS 
REPRESENTED  IN  FIG.  II 
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COMMENTS  AM)  RECOMMENDATIONS 


When  a  government  establishment  Is  put  In  an  isolated  area,  the 
regions  just  beyond  the  exits  may  soon  be  occupied  by  non-govemment 
employees  desiring  to  supply  goods  or  services  and  the  area  may  become 
a  site  for  housing  developments. 

If  possible,  a  test  center  involving  the  use  of  explosives  shotild 
be  so  located  that  the  prevailing  upper  vinds  will  blow  from  the  source 
of  blast  waves  toward  a  wide  stretch  of  Government  owned  lands  or  to  • 
the  open  sea. 

The  entrances  to  the  establishment  should  be  so  located  that  the 
prevailing  upper  winds  will  blow  from  them  toward  the  point  of  detonation. 

The  direction  of  the  winds  are  determined  by  the  position  of  the 
high  and  low  pressure  areas  in  the  atmosphere  and  with  a  little  experience 
in  this  procedure  a  study  of  the  daily  weather  maps  will  permit  a 
meteorologist  to  forecast  good  shooting  periods  several  days  in  advance. 
Since  the  conditions  in  the  atmosphere  usually  drift  toward  the  east, 
the  preceding  statement  will  be  of  no  help  to  stations  along  the 
Pacific  Coast  because  no  detailed  data  for  the  area  immediately  west  of 
the  coast  is  given  on  the  weather  maps. 

The  most  desirable  time  for  detonating  an  explosive  is  when  the 
center  of  a  high  pressure  area  is  at  the  point  of  detonation.  The 
temperature  of  the  air  around  the  test  site  will  then  be  found  to  decrease 
with  altitude  and  the  wind  will  be  at  a  minimum.  Itoder  such  conditions, 
the  blast  waves  are  refracted  upward  and  the  detonation  of  a  large 
charge  (1000  pounds  of  TNT)  would  be  barely  audible  at  a  distance  of  one 
or  two  miles. 

Generally  the  weather  conditions  change  slowly,  since  the  areas  of 
high  and  low  atmospheric  pressures  are  large  and  oh  the  average  move  only 
300  miles  per  day.  Therefore  the  conditions  indicated  by  the  meteor¬ 
ological  data  cein  be  expected  to  remain  sufficiently  constant  for  about 
8  hours.  The  firing  should  be  started  as  soon  as  possible  after  the 
meteorological  data  are  taken.  If  the  firing  is  to  extend  from  dayligdit 


51 


conditions  Into  the  night,  it  is  advisable  to  make  a  "metro"  run  in 
the  late  afternoon  to  detect  any  changes  that  are  occurring.  If  the 
weather  map  shows  a  "weather- front"  to  be  approaching  caution  in  the 
firing  is  advisable.  Firing  should  be  discontinued  during  the  passing 
of  the  "front"  and  a  metro  run  should  be  made  after  the  passage  of  the 
"front"  to  determine  if  firing  can  be  resumed. 

It  should  be  noted  that  the  worst  focusing  conditions  can  cause 
the  pressures  indicated  in  Figure  8  to  be  multiplied  by  100.  If  the., 
weight  of  the  explosive  being  detonated  is  so  small  that  the  pressure 
at  the  nearest  point  of  possible  damage  could  not  exceed  0.1  psl  even 
under  the  worst  conditions,  then  firing  can  proceed  without  a  metro 
run.  This  procedure  will  scmetlmes  cause  complaints  of  annoyance  but 
should  not  result  in  damage. 

RECAPITULATION 

The  following  summary  of  the  Standard  Operating  Procedure  will 
serve  as  a  convenient  guide; 

1.  The  meteorological  balloon  should  be  released  S  or  5  hours 
before  the  firing  program  Is  scheduled  to  begin. 

2.  The  wind  and  temperature  data  should  be  reduced  and  tabiilated 
as  in  Figure  10.  Data  should  be  recorded  at  intervals  of -500  feet  frcmi 
the  sijtrface  to  5000  feet  altitude  and  at  intervals  of  1000  feet  from 
5000  feet  to  12000  feet  altitude. 

5.  The  velocity  versus  altitude  shovO-d  be  plotted  for  each  azimuth 
of  interest  as  in  Figure  9  “the  velocity  gradients  should  be  noted, 

4.  The  combination  of  gradients  for  each  azimuth  should  be  compared 
to  the  set  of  32  in  Appendix  B  and  the  location  of  the  focus  (if  one 
is  indicated)  should  be  noted. 

(if  no  sample  velocity  versus  altitude  curve  in  Appendix  B  can  be 
found  to  match  the  one  just  plotted  then  the  analogue  computer  must  be 
used. ) 

.  5.  With  the  weight  of  explosives  to  be  detonated  and  the  dls-tance. 
to  the  focus  known,  the  pressure  to  be  expected  In  a  uniform  atmosphere 
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can  be  read  from  the  chart  on  Figure  8  of  the  text. 


6.  By  comparing  the  velocity  versus  altitude  curve  to  the  five 
basic  types  shown  in  Figure ■9»  the  multiplication  factor  for  the 
prevailing  conditions  can  be  determined. 


7.  If  the  pressiire  read  from  Figure  8,  multiplied  by  the 
multiplication  factor  determined  from  Figure  9j  Is  0-1  Psi  or  greater, 
firing  should  be  postponed.- 


PAUL  H.  LORRAIN 


WILLIAM  H.  TOWNSEiro  • 
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APFEpiX  A 


APPROXIMATIONS  IN  ANALOGUE 
SOLUTION  OF  RAY  PATBB 
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The  following,  ms  fximished  by  Mr.  DeBow  Owen  of  Piedmont  Division^ 
Sperry-Rand  Corporation,  Charlottsvllle,  Virginia. 

Since  the  Ray  Tracer  computing  system  Integrates  in  the  time  domain, 
y  must  be  expressed  as  a  function  of  time; 

Sine  (1)  ' 

dt  ds  dt  y  y  '  ' 

In  which  0  =  angle  between  the  ray  and  the  horizontal  and  V  Is  the  . 

y  y 

velocity  of  sound,  at  altitude  Y. 


dy 


^  (V  Sin  0  ) 
dt  '  y  y' 


d0 


1  P  P 

-  §  (Cos'"  8  -  Sln'=^  e  )  3^ 


^  dQ 

“5  (cos  28  )  ^ 


(2) 


d0  dV 

Since  Cos  0^  =  CV^,  5^  =  ^ 

and  substituting  in  (2) 


dv 

■J  =  -Vy  (cos  28^  ^ 


(3) 


In  the  Ray  Tracer  the  assumption  Is  made  that  the  horizontal  distance 
traveled  Is  equal  to  a  constant  average  velocity  times  the  travel  time. 


l.e.,  ==  =  V*' 


V  =  constant 
av 
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then  y  =  g  g  =  tan 


and  *y  =  V  Sec  0 

•'  a.y  y  dt 
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(6) 


and  again  since 


c  ■ 

Cob  0 


-V^ 


y  = 


[%y 


2 

Y  Cos  0 

y  y 


5 

dy 


(7) 


The  gradients  commonly  experienced  in  the  lower  atmosphere  are 
small  (O  to  .040  ft/sec/ft.)  and  persist  for  only  a  few  thousand  feet. 
For  this  reason,  the  rays  that  return  to  the  surface  will  traverse  a 
path,  the  direction  of  which  will  have  a  maximum  inclination  to  the 
horizontal  less  than  12  degrees..  Therefore  a  further  approximation  has 


been  made  that  Cos 


-v® 


«y  = 


Then  y  * 


av 


dV 

dy 


(8) 
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Equation  (8)  is  used  in  the  ray  tracer  in  lieu  of  the  more  precise 
equation  (2l)  of  the  text.  A  comparison  of  the  exact  solution  of  these 
two  equations  will  provide  an  estimate  of  the  error  caused  hy  the  ap¬ 
proximations.  When  the  initial  slope  of  the  ray  is  eleven  degrees,  the 
difference  hetween  the  true  and  the  approximate  maximum  height  reached 
hy  the  ray  is  1.6  per  cent.  The  difference  in  the  true  range  and  the 
approximate  range  is  1.2  per  cent.  These  differences  are  insignificant.. 
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TYPICAL  VERTICAL  VELOCITY  GRADIENTS 


59 


In  the  following  pages,  there  is  presented  a  set  of  gradients  which', 
are  typical  .of  conditions  found  at  Aberdeen  Proving  Ground  over  a  two 
year  period.  The  velocity  at  each  altitude  was  determined  from-  the 

I 

temperature  and  wind  velocity. 

¥ith  each  velocity  versus  altitude  curve  the  resulting  sound  ray 
paths,  as  computed  on  the  Sperry-Rand  Electronic  Ray  Plotter,  are'  • 
shown . 
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